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Abstract 

Blood  is  a complex  suspension  that  demonstrates  several  non-Newtonian 
rheological  characteristics  such  as  deformation-rate  dependency,  viscoelas- 
ticity and  yield  stress.  In  this  paper  we  outline  some  issues  related  to  the 
non-Newtonian  effects  in  blood  circulation  system  and  present  modeling  ap- 
proaches based  mostly  on  the  past  work  in  this  field. 
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1 Introduction 

Blood  is  a heterogeneous  multi-phase  mixture  of  solid  corpuscles  (red  blood  cells, 
white  blood  cells  and  platelets)  suspended  in  a liquid  plasma  which  is  an  aque- 
ous solution  of  proteins,  organic  molecules  and  minerals  (refer  to  Figure  1).  The 
rheological  characteristics  of  blood  are  determined  by  the  properties  of  these  com- 
ponents and  their  interaction  with  each  other  as  well  as  with  the  surrounding 
structures.  The  blood  rheology  is  also  affected  by  the  external  physical  conditions 
such  as  temperature;  however,  in  living  organisms  in  general,  and  in  large  mammals 
in  particular,  these  conditions  are  regulated  and  hence  they  are  subject  to  minor 
variations  that  cannot  affect  the  general  properties  significantly  [1].  Other  physical 
properties,  such  as  mass  density,  may  also  play  a role  in  determining  the  blood 
overall  rheological  conduct.  The  rheological  properties  of  blood  and  blood  vessels 
are  affected  by  the  body  intake  of  fluids,  nutrients  and  medication  [2-5],  although 
in  most  cases  the  effect  is  not  substantial  except  possibly  over  short  periods  of  time 
and  normally  does  not  have  lasting  consequences. 

The  viscosity  of  blood  is  determined  by  several  factors  such  as  the  viscosity 
of  plasma,  hematocrit  level  (refer  to  Figures  2 and  3),  blood  cell  distribution, 
and  the  mechanical  properties  of  blood  cells  [6-8].  The  blood  viscosity  is  also 
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Figure  1:  Simplified  illustration  of  the  main  blood  composition.  The  number  and 
size  of  cells  shown  in  the  figure  are  for  demonstration  purposes  and  not  supposed 
to  reflect  their  proportion  or  relative  size  in  the  blood. 


affected  by  the  applied  deformation  forces,  extensional  as  well  as  shearing,  and  the 
ambient  physical  conditions.  While  the  plasma  is  essentially  a Newtonian  fluid, 
the  blood  as  a whole  behaves  as  a non-Newtonian  fluid  showing  all  signs  of  non- 
Newtonian  rheology  which  includes  deformation  rate  dependency,  viscoelasticity, 
yield  stress  and  thixotropy.  Most  non-Newtonian  effects  originate  from  the  red 
blood  cells  due  to  their  high  concentration  and  distinguished  mechanical  properties 
such  as  elasticity  and  ability  to  aggregate  forming  three-dimensional  structures  at 
low  deformation  rates  [1,  9]. 

Deep  understanding  of  the  blood  rheology,  which  includes  its  non-Newtonian 
characteristics,  is  important  for  both  diagnosis  and  treatment.  Looking  into  the 
existing  biological  literature  and  comparing  to  the  non-biological  literature,  such 
as  earth  science  studies,  it  is  obvious  that  the  non-Newtouian  phenomena  in  the 
blood  circulation  has  not  been  given  sufficient  attention  in  the  biological  studies. 
One  reason  is  the  complexity  of  the  biological  systems  which  makes  the  considera- 
tion of  non-Newtonian  effects  in  blood  circulation  more  difficult  to  handle.  Hence, 
to  simplify  the  flow  models  and  their  computational  implementation,  blood  is  gen- 
erally assumed  Newtonian  and  these  effects  are  ignored.  The  obvious  difficulties 
in  observation,  experimentation  and  measurement  related  to  in  vivo  blood  flow 
add  another  complicating  factor.  Moreover,  apart  from  some  rare  and  extreme 
pathological  states,  the  non-Newtonian  effects  in  blood  flow  are  relatively  mild  as 
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Figure  2:  Dependency  of  human  blood  viscosity  on  hematocrit  percentage  concen- 
tration on  a linear-linear  plot  [6,  10]. 

compared  to  the  non-Newtonian  effects  exhibited  by  typical  polymeric  systems  for 
instance.  This  makes  the  approximation  of  blood  as  a Newtonian  fluid  an  accept- 
able assumption  and  not  far  from  reality  in  a significant  part  of  the  circulatory 
system  under  normal  conditions. 

Many  theoretical,  numerical,  experimental  and  clinical  studies  of  non-Newtonian 
effects  in  blood  circulation  have  been  conducted  in  the  last  few  decades.  However, 
there  is  no  general  approach  in  tackling  this  problem  in  a systematic  way  based  on 
a unified  vision.  Almost  all  the  past  studies  focus  on  individual  problems  and  deal 
with  the  existing  non-Newtonian  phenomena  within  a limited  local  context.  The 
current  study,  which  is  basically  a brief  overview  of  this  subject,  is  trying  to  deal 
with  the  non-Newtonian  blood  rheology  in  general  as  applied  to  all  levels  of  the 
circulation  system. 

2 Non-Newtonian  Characteristics 

As  indicated  already,  blood  is  a complex  non-Newtonian  fluid  showing  various  signs 
of  non-Newtonian  rheology  such  as  shear  thinning  [11],  yield  stress  [7,  12-16]  and 
viscoelasticity  [1,  13,  14,  17-20].  The  blood  is  also  characterized  by  a distinctive 
thixotropic  behavior  [9,  21]  revealed  by  the  appearance  of  hysteresis  loops  during 
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shearing  cycles  [21-23].  These  non-Newtonian  properties  do  not  affect  the  flow 
patterns  inside  the  flow  paths  and  the  fluid  transportation  only  but  they  also  affect 
the  mechanical  stress  on  the  blood  vessel  walls  and  the  surrounding  tissues  [24] 
especially  in  cases  of  irregular  lumen  geometry  like  stenosed  arteries  [25-27].  The 
mechanical  stress  on  the  vessel  wall  and  tissue  is  not  only  important  for  its  direct 
mechanical  impact,  especially  when  sustained  over  long  periods  of  time,  but  it  can 
also  contribute  to  the  commencement  and  advancement  of  long  term  lesions  such  as 
forming  sediments  inside  the  vessel  wall  [28].  The  non-Newtonian  properties,  like 
viscoelasticity,  have  also  an  impact  on  other  transport  phenomena  such  as  pulse 
wave  propagation  in  arteries  [18]. 

Non-Newtonian  effects  in  general  are  dependent  on  the  magnitude  of  deforma- 
tion rates  and  hence  they  can  exist  or  be  enhanced  at  certain  flow  regimes  such  as 
low  shear  rates  [20,  29].  The  non-Newtonian  effects  are  also  influenced  by  the  type 
of  deformation,  being  shear  or  elongation  [16].  The  impact  of  the  non-Newtonian 
effects  can  be  amplified  by  a number  of  factors  such  as  pathological  blood  rheology 
and  flow  in  stenosed  vessels  and  stents  [30-33].  An  interesting  finding  of  one  study 
[34]  is  that  although  flow  resistance  and  wall  shear  stress  increases  as  the  size  of 
stenosis  increases  for  a given  rheological  model,  the  non-Newtonian  nature  of  blood 
acts  as  a regulating  factor  to  reduce  the  resistance  and  stress  and  hence  contribute 
to  the  body  protection.  In  this  context,  shear  thinning  seems  to  have  the  most 
significant  role  in  facilitating  blood  flow  through  stenotic  vessels. 

Blood  is  a predominantly  shear  thinning  fluid  (refer  to  Figure  3),  especially 
under  steady  flow  conditions,  and  this  property  has  the  most  important  non- 
Newtonian  impact  [7,  20,  24,  35].  Shear  thinning  is  not  a transient  characteris- 
tics; moreover  it  is  demonstrated  at  most  biological  flow  rates  [35]  although  it  is 
more  pronounced  at  low  deformation  regimes.  Shear  thinning  rheology  arises  from 
disaggregation  of  the  red  blood  cells  with  increasing  shear  rate  [36].  This  same  rea- 
son is  behind  the  observed  thixotropic  blood  behavior  as  shearing  forces  steadily 
disrupt  the  structured  aggregation  of  blood  cells  with  growing  deformation  time. 
The  origin  of  other  non-Newtonian  effects  can  also  be  traced  back  to  the  blood 
microstructure  as  will  be  discussed  next. 

The  viscoelastic  nature  of  blood  basically  arises  from  its  corpuscular  microstruc- 
ture.  Viscoelastic  properties  originate  from  the  red  blood  cells,  which  are  distin- 
guished by  their  pronounced  elastic  deformability  [36,  39]  associated  with  the  abil- 
ity to  aggregate  forming  three-dimensional  structures  known  as  rouleaux  [30].  The 
aggregation  is  mostly  demonstrated  at  low  shear  rates  and  hence  non-Newtonian 
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Figure  3:  Human  blood  viscosity  as  a function  of  shear  rate  for  a range  of  hemat- 
ocrit concentrations  on  a log-log  plot  [8,  37,  38]. 

behavior  in  general  and  viscoelasticity  in  particular  are  more  pronounced  at  these 
regimes  of  low  deformation  [1,  6,  40].  The  viscoelastic  effects  are  magnified,  if  not 
activated,  by  the  pulsatile  nature  of  blood  flow.  The  viscoelastic  effects  in  blood 
circulation  should  not  be  limited  to  the  viscoelastic  properties  of  the  blood  itself 
but  also  to  the  viscoelastic  (or  elastic  depending  on  the  adopted  model)  properties 
of  the  blood  vessels  and  the  porous  tissue  through  which  the  blood  is  transported 
[41],  This  can  be  justified  by  the  fact  that  all  these  effects  are  manifested  by  the 
blood  circulation  and  hence  they  participate,  affect  and  affected  by  the  circulation 
process.  An  interaction  between  the  viscoelastic  behavior  of  blood  with  that  of  the 
vessel  wall  and  porous  tissue  is  unavoidable  consequence. 

Blood  also  demonstrates  yield  stress  although  there  is  a controversy  about  this 
issue  [1,  25].  Yield  stress  arises  from  the  aggregation  of  red  blood  cells  at  low  shear 
rates  to  form  the  above-mentioned  three-dimensional  micro-structures  (rouleaux) 
that  resist  the  flow  [1,  7,  42],  Studies  have  indicated  that  yield  stress  is  posi- 
tively correlated  to  the  concentration  of  fibrinogen  protein  in  blood  plasma  and 
to  the  hematocrit  level  [17,  25,  39,  43].  An  illustrative  plot  of  the  dependence  of 
yield  stress  on  hematocrit  level  is  shown  in  Figure  4.  Other  factors,  such  as  the 
concentration  of  minerals,  should  also  have  a contribution.  Many  of  the  blood 
rheological  characteristics  in  general,  and  non-Newtonian  in  particular,  are  also 
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controlled  or  influenced  by  the  fibrinogen  level  [42],  The  yield  stress  characteristic 
of  blood  seems  to  vanish  or  become  negligible  when  hematocrit  level  falls  below  a 
critical  value  [17].  Yield  stress  contributes  to  the  blood  clotting  following  injuries 
and  subsequent  healing,  and  may  also  contribute  to  the  formation  of  blood  clots 
(thrombosis)  and  vessel  blockage  in  some  pathological  cases  such  as  strokes.  The 
value  of  yield  stress,  as  reported  in  a number  of  clinical  and  experimental  studies, 
seems  to  indicate  that  it  is  not  significant  and  hence  it  has  no  tangible  effect  on 
the  flow  profile  (and  hence  flow  rate)  at  the  biological  flow  ranges  in  large  and 
medium  size  blood  vessels  [25].  However,  it  should  have  more  significant  impact 
in  the  minute  capillaries  and  some  porous  structures  where  flow  at  very  low  shear 
rates  occurs.  The  magnitude  of  yield  stress  and  its  effect  could  be  aggravated  by 
certain  diseased  states  related  to  the  rheology  of  blood,  like  polycythemia  vera,  or 
the  structure  of  blood  vessels  such  as  stenoses. 


Figure  4:  Dependence  of  yield  stress  of  human  blood  on  hematocrit  level  on  a log- 
linear  plot  according  to  some  studies.  Other  types  of  correlation  have  also  been 
obtained  in  other  studies  [8,  14,  44,  45]. 

As  a shear  thinning  fluid,  blood  is  also  characterized  by  a thixotropic  behavior, 
which  is  confirmed  experimentally  by  a number  of  studies  [9],  due  to  the  intimate 
relation  between  these  two  non-Newtonian  properties  [23].  This  may  also  explain  a 
possible  controversy  about  the  thixotropic  nature  of  blood  [46]  as  the  thixotropic- 
like  behavior  may  be  explained  by  other  non-Newtonian  characteristics  of  blood. 
Despite  the  fact  that  thixotropy  is  a transient  property,  due  to  the  pulsative  nature 
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of  the  blood  flow  the  thixotropic  effects  may  have  long  term  impact  on  the  blood 
circulation.  This  equally  applies  to  the  time-dependent  effects  of  viscoelasticity. 
Thixotropy  is  more  pronounced  at  low  shear  rates  with  a long  time  scale.  The 
effect,  however,  seem  to  have  a less  important  role  in  blood  flow  than  other  non- 
Newtonian  effects  such  as  shear  thinning  [35],  and  this  could  partly  explain  the 
limited  amount  of  studies  dedicated  to  this  property.  The  thixotropic  behavior  of 
blood  is  very  sensitive  to  the  blood  composition  and  hence  it  can  demonstrate  big 
variations  between  different  individuals  and  under  different  biological  conditions 
[9], 

It  should  be  remarked  that  time  dependent  effects  in  general,  whether  thixotropic 
or  viscoelastic  in  nature  or  of  any  other  type,  should  be  expected  in  the  flow  of 
blood  due  to  a number  of  reasons.  One  reason  is  the  pulsatility  of  blood  flow  and 
the  rapid  change  in  the  deformation  conditions  during  the  systolic-diastolic  cardiac 
cycle.  Another  reason  is  the  rapid  change  in  the  shear  magnitude  between  one  part 
of  the  system  to  another  part,  i.e.  different  shear  rates  between  the  arteries,  capil- 
laries, porous  tissue  and  venous  part  [47].  A third  reason  is  the  irregular  shape,  such 
as  bends  and  converging-diverging  formations,  of  the  blood  flow  conduits  [48,  49] 
which  activates  or  accentuates  time-dependent  effects.  A fourth  reason  is  the  dif- 
ference in  the  deformation  rates  between  the  ventricular  systole  and  diastole.  The 
last  reason  may  explain  the  indication  of  one  study  [7]  that  the  non-Newtonian 
effects  are  more  important  at  diastole  than  systole  since  the  shear  rates  during 
diastole  are  expected  to  be  lower  than  those  at  systole. 

Another  remark  is  that  most  of  the  reported  non-Newtonian  rheological  pa- 
rameters, as  well  as  many  other  physical  properties  of  blood,  are  obtained  from  in 
vitro  measurements  and  hence  they  are  subject  to  significant  errors  as  an  indicator 
to  the  in  vivo  values  due  to  the  difference  in  ambient  conditions  as  well  as  the 
experimental  requirements  and  procedures,  such  as  using  additives  to  preserve  and 
fluidize  the  blood  samples,  that  can  introduce  significant  variations  on  the  blood 
properties.  Moreover,  the  reported  values  could  be  highly  dependent  on  the  mea- 
surement method  [8].  The  differences  between  the  individual  subjects  and  their 
conditions  like  dietary  intake  prior  to  measurement  [3],  most  of  which  are  difficult 
to  control  or  quantify,  should  add  to  the  uncertainties  and  fluctuations.  Hence, 
most  of  these  values  should  be  considered  with  caution  especially  when  used  for  in 
vivo  and  patient-specific  modeling  and  investigation. 
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3 Modeling  non-Newtonian  Effects 


Blood  is  a complex  non-Newtonian  fluid  and  hence  reliable  modeling  of  blood  flow 
in  the  circulation  system  should  take  into  account  its  non-Newtonian  character- 
istics. Several  non-Newtonian  rheological  models  have  been  used  to  describe  the 
blood  rheology;  these  models  include  Carreau-Yasuda  [7,  11,  24,  26,  50-58],  Casson 
[4,  7,  8,  14,  17,  25,  29,  33,  34,  37,  43,  54,  58,  59],  power  law  [1,  7,  20,  26,  34,  58,  60], 
Cross  [52-54,  58,  61],  Herschel-Bulkley  [8,  27,  56,  62,  63],  Oldroyd-B  [1,  36,  40,  64], 
Quemada  [4,  54,  58,  65],  Yeleswarapu  [55,  57,  66],  Bingham  [25],  Eyring- Powell 
[67],  and  Ree- Eyring  [4],  The  constitutive  equations  of  these  rheological  models 
are  given  in  Table  1.  Other  less  known  fluid  models  [31]  have  also  been  used  to 
describe  the  rheology  of  blood.  A quick  inspection  of  the  blood  literature  reveals 
that  the  most  popular  models  in  non-Newtonian  hemorheologic  and  hemodynamic 
modeling  are  Carreau-Yasuda  and  Casson. 

Blood  is  also  modeled  as  a Newtonian  fluid  [1,  26,  29,  60]  which  is  a good 
approximation  in  many  circumstances  such  as  the  flow  in  large  vessels  at  medium 
and  high  shear  rates  under  non-pathological  conditions  [59].  As  there  is  no  sudden 
transition  from  non-Newtonian  to  Newtonian  flow  as  a function  of  shear  rate,  there 
is  no  sharply-defined  critical  limit  for  such  a transition  [66]  and  hence  this  remains 
a matter  of  choice  which  depends  on  a number  of  objective  and  subjective  factors. 
However,  there  seems  to  be  a general  consensus  that  the  shear  rate  range  for  which 
non-Newtonian  effects  are  considered  significant  is  < 100  s-1  [6,  25,  42,  58,  60]; 
above  this  limit  the  blood  is  generally  treated  as  a Newtonian  liquid. 

No  single  model,  Newtonian  or  non-Newtonian,  can  capture  all  the  features 
of  the  blood  complexities  [66]  and  hence  different  models  are  used  to  represent 
different  characteristics  of  the  blood  rheology.  These  models,  whether  Newtonian 
or  non-Newtonian,  obviously  have  significant  differences  and  hence  they  can  pro- 
duce very  different  results  [55,  57,  61].  The  results  also  differ  significantly  between 
Newtonian  and  non-Newtonian  models  in  most  cases  [53,  58].  The  non-Newtonian 
models  vary  in  their  complexity  and  ability  to  capture  different  physical  phenom- 
ena. 

Diverse  methods  have  been  used  in  modeling  and  simulating  non-Newtonian 
effects  in  blood  rheology;  these  include  analytical  [14,  27,  63],  stochastic  [33,  59,  68- 
70],  and  numerical  mesh  methods;  such  as  finite  element  [24,  41,  51,  61,  71],  finite 
difference  [25,  31,  32,  40,  57,  64],  finite  volume  [1,  7,  29,  52,  53,  57,  58],  and  spectral 
collocation  methods  [40,  64], 

As  indicated  previously,  most,  if  not  all,  non-Newtonian  characteristics  arise 
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Table  1:  The  non-Newtonian  fluid  models  that  are  commonly  used  to  describe  blood 
rheology.  The  meanings  of  symbols  are  given  in  Nomenclature  § 6.  The  symbols 
that  define  fluid  characteristic  properties,  such  as  A,  are  generically  used  and  hence 
the  same  symbol  may  represent  different  physical  attributes  in  different  models. 
Some  of  these  models  may  have  more  than  one  form;  the  one  used  in  this  table  is 
the  widespread  of  its  variants.  The  last  column  represents  the  frequently  obtained 
non-Newtonian  properties  from  these  models  in  the  context  of  blood  modeling 
although  other  properties  may  also  be  derived  and  employed  in  modeling  other 
materials. 


Model 

Equation 

Non-Newtonian  Properties 

Carreau-Yasuda 

II  — II  4-  Mo- M°o 

P'  P'00  1 1 — n 

[1+(A7)°1  — 

shear  thinning 

Casson 

= ( kj Y + To*2 

yield  stress 

Power  law 

r = 

shear  thinning 

Cross 

ii  — ij  4-  w_^°° 

h1  P'00  1 l+X^m 

shear  thinning 

Herschel-Bulkley 

T = + T0 

shear  thinning,  yield  stress 

Oldroyd-B 

T + Air  = /i()  ^7  + A 27^ 

viscoelasticity 

Quemada 

( 1 ko+kooW^/ic  \ 

"“p1  2(1+VVPJ 

shear  thinning 

Yeleswarapu 

H ~ Hoo+  Po  ACo)  1+A^ 

shear  thinning 

Bingham 

r = k'j  + t0 

yield  stress 

Eyring-Powell 

, (M0-Moo)sinh“1(A7) 

p-  P'00  1 ^ 

shear  thinning 

Ree-Eyring 

t = tc  sinh-1  (^f1  j 

shear  thinning 

from  the  blood  microstructure  and  particularly  the  concentration,  distribution  and 
mechanical  properties  of  the  red  blood  cells.  For  example,  the  viscoelastic  prop- 
erties of  blood  originate  from  the  mechanical  properties  of  the  suspended  cells 
and  their  capability  of  elastic  deformation  and  structural  formation,  while  the 
thixotropic  properties  arise  from  steady  disaggregation  of  blood  cells  over  pro- 
longed shearing  time.  Interestingly,  the  majority  of  the  rheological  models  used  to 
describe  the  blood  rheology  are  bulk  phenomenological  models  of  empirical  nature 
with  little  consideration,  if  any,  to  its  highly  influential  micro-structure.  Hence, 
more  structurally-based  models,  such  as  those  based  on  molecular  dynamics,  are 
required  to  improve  the  description  and  modeling  of  the  blood  rheological  behavior. 
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4 Non-Newtonian  Effects  in  Circulation  Subsys- 
tems 

Because  the  impact  of  the  non-Newtonian  effects  is  highly  dependent  on  the  shape 
and  size  of  the  flow  conduits,  different  non-Newtonian  rheological  behavior,  and 
hence  different  flow  modeling  approaches,  should  apply  to  the  different  parts  of  the 
circulatory  system.  Different  approaches  are  also  required  because  of  the  difference 
in  the  nature  of  the  blood  transportation  processes  in  these  parts,  such  as  large 
scale  bulk  flow  in  the  large  vessels  as  opposite  to  perfusion  or  diffusion  in  the  porous 
tissue. 

We  can  identify  three  types  of  circulatory  subsystems  in  which  non-Newtonian 
effects  should  be  analyzed  and  modeled  differently:  large  blood  vessels  which  mainly 
apply  to  arteries  and  veins,  small  blood  vessels  which  broadly  include  capillaries 
and  possibly  arterioles,  and  porous  tissue  such  as  the  myocardium  and  muscles 
in  general.  These  three  subsystems  are  graphically  illustrated  in  Figure  5.  The 
distinction  between  large  and  small  vessels  is  not  clear  cut  as  it  depends  on  the  na- 
ture of  the  flow  phenomenon  under  consideration  and  the  associated  circumstances. 
However,  the  distinctive  feature  that  should  be  used  as  a criterion  to  differentiate 
between  these  two  categories  of  blood  vessels  in  this  context  is  the  validity  of  the 
continuum  approximation  as  applied  to  the  blood  where  in  the  large  vessels  this 
approximation  is  strictly  held  true  while  in  the  small  vessels  it  approaches  its  limits 
as  some  non-continuum  phenomena  start  to  appear. 

In  the  following  subsections  we  outline  general  strategies  for  modeling  non- 
Newtonian  effects  in  the  circulation  subsystems. 

4.1  Large  Vessels 

In  large  vessels,  which  include  large  cavities  such  as  the  ventricles  and  atria  inside 
the  myocardium  as  well  as  the  large  arteries  and  veins,  the  blood  essentially  behaves 
as  a Newtonian  fluid.  One  reason  is  that  the  blood  in  such  large  lumens  and 
cavities  is  normally  exposed  to  a relatively  high  shear  rates  and  hence  the  non- 
Newtonian  effects  which  are  basically  induced  at  low  shear  rates  die  out  [36].  Also 
at  this  large  scale  the  blood  appears  as  a homogeneous  continuum  medium  with 
diminishing  effect  of  blood  cell  aggregation  [50] . The  interaction  between  the  blood 
cells  with  their  pronounced  elastic  properties  is  also  minimal  at  this  scale.  However, 
in  some  pathological  situations,  non-Newtonian  effects  are  important  even  in  the 
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(a)  Large  Vessels 


(b)  Small  Vessels 


(c)  Porous  Tissue 


Figure  5:  The  three  types  of  circulatory  subsystems:  (a)  large  vessels,  (b)  small  ves- 
sels, and  (c)  porous  tissue.  The  three  sub-figures  are  obviously  drawn  on  different 
scales. 
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big  cavities  and  large  vessels  and  therefore  they  should  be  considered  in  the  flow 
model.  This  may  also  be  true  in  some  non-pathological  situations  in  which  the 
non-Newtonian  effects  can  be  critical  to  the  observed  phenomena. 

ft  should  be  remarked  that  the  non-Newtonian  effects  in  the  venous  part  of  the 
circulatory  system  should  be  more  important  than  the  arterial  part  due  to  the  lower 
deformation  rates  in  the  former  than  the  latter  as  this  seems  to  be  an  accepted  fact 
in  the  hemodynamic  studies,  ffowever,  we  did  not  find  a proper  discussion  about 
this  issue  in  the  available  literature.  The  difference  in  the  blood  composition  in 
these  two  parts;  due  for  example  to  the  difference  in  concentration  of  substances 
like  nutrients,  oxygen  and  metabolic  wastes;  should  also  affect  the  non-Newtonian 
rheology  in  these  two  subsystems  and  introduce  more  complications  in  modeling 
blood  flow  especially  in  large  vessels.  Like  the  previous  issue,  we  did  not  find  an 
explicit  discussion  to  this  issue  in  the  available  literature.  Another  remark  is  that 
certain  parts  of  the  large  vessels  network  can  contain  spots  of  low  shear  rates  such 
as  bends  and  bifurcation  junctions  [7,  25]  and  hence  non-Newtonian  effects  in  large 
vessels  could  be  significant  in  some  cases  where  these  spots  play  an  exceptionally 
important  role  in  the  blood  flow  due  to  a diseased  case  for  instance. 

Several  mathematical  and  computational  models  have  been  used  to  describe 
the  flow  of  blood  in  large  individual  vessels.  These  models  include  the  elastic  one- 
dimensional Navier-Stokes  and  the  rigid  Hagen-Poiseuille  for  Newtonian  fluids,  as 
well  as  many  other  non-Newtonian  rheological  models  such  as  Cross  and  Carreau- 
Yasuda,  as  discussed  previously  is  section  3.  The  characteristics  of  blood  flow  in 
large  single  vessels  are  obtained  from  these  mathematical  models  either  analytically 
or  numerically,  e.g.  through  the  use  of  finite  element  or  finite  difference  techniques. 
Most  of  the  employed  non-Newtonian  fluid  models  are  generalized  Newtonian  mod- 
els and  hence  they  do  not  account  for  history-dependent  elastic  or  thixotropic  ef- 
fects. Also,  the  analytical  non-Newtonian  models  generally  apply  to  rigid  tubes 
only  although  there  are  some  attempts  in  this  context  to  extend  Poiscuillc  flow  to 
elastic  vessels  with  non-Newtonian  rheology  [63].  Numerical  methods  may  also  be 
used  to  extend  the  non-Newtonian  models  to  elastic  vessels. 

With  regards  to  the  flow  in  vascular  networks  of  large  vessels,  the  main  models 
used  to  describe  and  simulate  blood  flow  are  the  one-dimensional  Navier-Stokes 
finite  element  model  for  elastic  networks  [72,  73]  and  the  Hagen-Poiseuille  model 
for  rigid  networks  [74],  Both  of  these  models  are  Newtonian  although  the  second 
may  be  extended  to  Poiseuillc-like  non-Newtonian  flow  through  the  inclusion  of 
time-independent  non-Newtonian  effects  using  a vessel-dependent  non-Newtonian 
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effective  viscosity  which  is  computed  and  updated  iteratively  to  reach  a consistent 
flow  solution  over  the  whole  network  [16,  75-77].  Based  on  a non-thorough  inspec- 
tion of  the  available  literature,  there  seems  to  be  no  extension  to  the  traditional 
one- dimensional  Navier-Stokes  distensible  network  model  [72,  78,  79]  to  incorporate 
non-Newtonian  effects.  In  fact  we  did  not  find  serious  attempts  in  the  available  lit- 
erature to  extend  the  Navier-Stokes  equation  in  general  (whether  one-dimensional 
or  multi-dimensional,  for  rigid  or  elastic  conduits)  to  account  for  non-Newtonian 
effects,  although  there  are  attempts  to  incorporate  non-Newtonian  effects  numer- 
ically into  Navier-Stokes  flow  models  [80,  81].  The  Navier-Stokes  equation  with 
its  nonlinearity  is  sufficiently  complex  to  be  solved  for  Newtonian  flow  in  most 
cases  let  alone  with  the  added  complexities  and  nonlinearities  introduced  by  the 
non-Newtonian  rheology. 

With  regards  to  the  one-dimensional  Navier-Stokes  distensible  model  for  single 
vessels  and  networks,  we  propose  two  possible  general  approaches  for  extending 
this  model  to  account  for  non-Newtonian  effects.  One  approach  is  to  accommo- 
date these  effects  in  the  fluid  viscosity  as  parameterized  by  the  viscosity  friction 
coefficient  k.  A second  possible  approach  is  to  incorporate  these  effects  in  the 
flow  profile  as  described  by  the  momentum  flux  correction  factor  a in  the  one- 
dimensional model  [72,  73].  For  the  network,  the  second  approach  is  based  on 
defining  a vessel- dependent  a field  over  the  whole  network.  A similar  viscosity 
field  may  also  be  required  for  the  first  approach  as  well. 

4.2  Small  Vessels 

Non-Newtonian  effects  are  generally  more  pronounced  in  small  flow  ducts,  such 
as  capillaries,  than  in  large  ducts  like  arteries  due  to  several  reasons  such  as  the 
deterioration  of  the  continuum  assumption  at  small  scales  especially  for  complex 
dispersed  systems  like  blood.  In  such  vessels  the  continuum  approximation  reaches 
its  limit  and  the  effect  of  blood  cell  aggregation  with  their  interaction  with  the 
vessel  wall  becomes  pronounced.  This  activates  the  non-Newtonian  rheological 
flow  modes  such  as  the  induction  of  elastic  effects  which  are  associated  with  the 
elastic  properties  of  the  red  blood  cells  and  their  structural  formation.  Also,  the 
non-Newtonian  effects  of  blood  are  more  prominent  at  low  shear  rates  which  are  the 
dominant  flow  regimes  in  the  small  vessels  [54],  Hence  non-Newtonian  rheological 
effects  should  be  considered  in  modeling,  simulating  and  analyzing  the  flow  of  blood 
in  small  vessels. 
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4.3  Porous  Tissue 


The  commonly  used  approach  in  modeling  blood  perfusion  in  living  tissue  is  to 
treat  the  tissue  as  a spongy  porous  medium  and  employ  Darcy  law  [23,  71,  82- 
85]  which  correlates  the  volumetric  flow  rate  to  the  pressure  gradient.  There  are 
several  limitations  in  the  Darcy  flow  model  in  general  and  in  its  application  to  the 
blood  flow  through  biological  tissue  in  particular,  and  hence  remedies  have  been 
proposed  and  used  to  improve  the  model.  Since  this  law  is  originally  developed  for 
the  flow  through  rigid  porous  media,  modified  versions  are  normally  used  to  account 
for  elasticity  as  required  for  modeling  biological  tissue  [71,  86].  Other  limitations 
include  neglecting  edge  effects  and  the  restriction  imposed  by  the  laminar  low 
velocity  assumption  on  which  the  Darcy  flow  is  based.  The  former  may  be  overcome 
by  employing  the  boundary  term  in  the  Brinkman  equation  while  the  latter  can  be 
eliminated  through  the  use  of  the  Forchheimer  model  which  incorporates  a high- 
velocity  inertial  term  [83]. 

Since  Darcy  law  in  its  original  formulation  is  based  on  the  Newtonian  flow  as- 
sumptions, non-Newtonian  rheology  is  generally  ignored  in  the  modeling  of  blood 
perfusion  through  porous  tissue.  There  have  been  several  extensions  and  modifica- 
tions to  the  Darcy  law  to  include  non-Newtonian  effects  in  the  flow  of  fluids  in  gen- 
eral, and  polymers  in  particular,  through  rigid  non-biological  porous  media.  These 
attempts  include,  for  example,  viscoelastic  models  [87-90],  Herschel-Bulkley  [91], 
power  law  [92,  93],  Blake-Kozeny-Carman  [94],  as  well  as  other  non-Newtonian  pro- 
totypes [95,  96].  Pore  scale  network  modeling  has  also  been  used  to  accommodate 
various  non-Newtonian  effects;  such  as  shear  thinning,  yield  stress  and  viscoelas- 
ticity; in  the  flow  of  polymers  through  rigid  porous  media  [15,  19,  75-77,  97-100]. 
Similarly,  other  computational  techniques,  such  as  stochastic  lattice  Boltzmann 
[101,  102],  have  been  tried  to  simulate  and  investigate  the  non-Newtonian  effects 
of  the  flow  through  rigid  porous  media  in  non-biological  studies.  However,  it  seems 
there  is  hardly  any  work  on  modeling  the  non-Newtonian  effects  in  the  blood  flow 
through  living  tissues  by  incorporating  these  effects  into  the  distensible  Darcy  flow 
model. 

To  conclude,  non-Newtonian  effects  in  the  blood  perfusion  through  porous  tis- 
sue are  not  negligible  in  general  due  to  the  existence  of  fluid  shearing  and  exten- 
sional  forces  which  activate  non-Newtonian  rheology.  As  the  deformation  rates 
in  this  type  of  transportation  is  normally  low,  and  considering  the  small  size  and 
tortuous  converging-diverging  shape  of  the  porous  space  inside  which  the  blood 
perfuses,  non-Newtonian  rheological  effects  are  expected  to  be  significant.  Non- 
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Newtonian  rheological  effects  associated  with  other  fluid  transport  phenomena, 
such  as  diffusion,  could  be  negligible  in  such  porous  space  due  to  the  absence  of 
shearing  and  extensional  forces  as  a result  of  lack  of  large  scale  fluid  bulk  move- 
ment in  these  micro-  and  nano-scale  phenomena.  However,  the  causes  underlying 
the  non-Newtonian  rheology  should  affect  these  transport  phenomena  as  well,  al- 
though more  serious  investigations  are  required  to  reach  any  definite  conclusion 
about  these  issues.  The  existing  literature  is,  unfortunately,  limited  in  this  scope 
[103,  104], 

5 Conclusions 

Blood  is  essentially  a non-Newtonian  suspension.  Its  complex  rheological  non- 
Newtonian  behavior  is  largely  influenced  by  its  microstructure  which,  through  the 
essentially  viscous  water-based  Newtonian  plasma  combined  with  the  effect  of  ag- 
gregation, deformation  and  orientation  of  the  suspended  blood  cells  with  their  dis- 
tinguished elastic  and  three-dimensional  structural  formation  properties,  can  show 
diverse  non-Newtonian  effects  at  various  shear  rate  regimes  and  through  different 
flow  conduit  structures  although  these  effects  are  more  pronounced  at  certain  flow 
regimes  and  in  particular  structures. 

Blood  rheological  properties,  and  its  mechanical  characteristics  in  general,  are 
affected  by  several  factors  such  as  the  type  and  magnitude  of  deformation  rate, 
hematocrit  level  and  protein  concentration.  Because  blood  is  a suspension,  its 
properties  can  be  strongly  influenced  by  the  shape  and  size  of  its  flow  conduits. 
The  non-Newtonian  effects  of  blood  can  be  accentuated  by  certain  pathological 
conditions  such  as  hypertension  and  myocardial  infarction. 

Apart  from  some  extreme  diseased  states  where  the  non-Newtonian  effects 
play  an  exceptionally  important  role  in  the  fluid  transport  phenomena,  the  non- 
Newtonian  effects  are  generally  mild  in  the  bulk  flow  of  blood  in  large  vessels.  More 
important  influence  of  non-Newtonian  rheology  occurs  in  the  flow  of  blood  through 
small  vessels  and  in  the  blood  perfusion  through  porous  tissue.  Other  transport 
phenomena,  like  diffusion,  which  do  not  involve  bulk  flow  associated  with  deforma- 
tion forces  of  shearing  or  extensional  type  that  activate  non-Newtonian  rheology, 
should  not  be  affected  directly  by  the  non-Newtonian  characteristics  although  the 
physical  causes  at  the  root  of  the  non-Newtonian  rheology  should  have  an  im- 
pact on  these  processes.  More  fundamental  studies  are  required  to  reach  specific 
conclusions  about  these  issues. 
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Non-Newtonian  rheology,  such  as  viscoelasticity,  of  the  blood  vessel  walls  and 
the  spongy  porous  tissue  should  also  be  considered  as  a contributor  to  the  over- 
all non-Newtonian  behavior  of  blood  circulation  as  these  effects  are  both  non- 
Newtonian  and  circulatory  in  nature  like  the  ones  demonstrated  by  the  blood  it- 
self. The  effect  of  fluid-structure  interaction  should  also  be  included  in  analyzing, 
modeling  and  simulating  of  non-Newtonian  effects  in  the  blood  circulation  as  it 
plays  an  important  hemodynamic  and  hemorheologic  role. 
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6 Nomenclature 

7 shear  rate 

7C  characteristic  shear  rate 

7 rate  of  strain  tensor 

A characteristic  time  constant 

Ai  relaxation  time 

A2  retardation  time 

/z  fluid  viscosity 

Up  plasma  viscosity 

/io  zero-shear-rate  viscosity 

/Zoo  infinite-shear-rate  viscosity 

r shear  stress 

r stress  tensor 

rc  characteristic  shear  stress 

t0  yield-stress 

(j)  volume  concentration 

a Carreau-Yasuda  index 

k consistency  coefficient 

k0  maximum  volume  fraction  for  zero  shear  rate 
koo  maximum  volume  fraction  for  infinite  shear  rate 
m Cross  model  index 

n power  law  index 

upper  convected  time  derivative 
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